%% test TR %%
Td = 38;
T0 = 35.4;
t = 0:0.1:Td;
Tr = zeros(1, 10*Td+1);
for i = 1:(10*Td+1)
    Tr(i) = TR(Td, T0, t(i));
end
plot(t, Tr, 'LineWidth',2);